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We study the spin-parity distribution P(J 7r ,E) of 156 Gd excited states above the neutron separation 
energy that are expected to be populated via the neutron pickup reaction 157 Gd( 3 He, 4 He) 156 Gd. In 
general, modeling of the spin-parity distribution is important for the applicability of the surrogate 
reaction technique as a method of deducing reaction cross sections. We model excited states in 156 Gd 
as rotational states built on intrinsic states consisting of a hole in the core where the hole represents 
neutron removal form a deformed single particle state. The reaction cross section to each excited 
state is calculated using standard reaction code that uses spherical reaction form- factor input. The 
spectroscopic factor associated with each form-factor is the expansion coefficient of the deformed 
neutron state in a spherical Sturmian basis consisting of the spherical reaction form-factors. 

PACS numbers: 21.10.Jx, 24.50.+g,24.10.Eq,25.55.Ci 



I. INTRODUCTION AND MOTIVATION 

Most light nuclei up to C and O are produced via nu- 
clear reactions within the stars. Heavier nuclei near the 
valley of stability are produced via slow neutron capture 
process (s process). To account for the observed abun- 
dance of the elements, however, one needs to also con- 
sider other processes such as proton capture process (p 
process) and rapid neutron capture process (r process). 
The r process is essential in the production of the heav- 
iest neutron rich elements since this process generates 
nuclei far from the valley of stability that decay back to- 
wards the valley. Successively heavier neutron rich nuclei 
are produced through violent processes such as supernova 
explosions. Such nuclei usually decay quickly. As a con- 
sequence, the understanding of the r-process contribu- 
tion to the observed abundance of the elements requires 
knowledge of the neutron induced reactions on unstable 
nuclei. Unfortunately, measuring these cross sections in a 
laboratory environment is a very difficult, if not impossi- 
ble, task because of the technical and practical problems 
associated with the use of an unstable nuclei. 

There have been various proposals for circumventing 
the problems associated with the use of unstable nuclei 
and yet to gain information about the desired nuclear 
reaction. One such approach is the Surrogate Method 
shown schematically in Fig. [T] 

In this paper we study the reaction 3 He+ 157 Gd — > 
4 He+ 156 Gd* and model the formation probability of var- 
ious excited states of the 156 Gd system within a direct 
reaction framework. This reaction is a testing ground 
of the surrogate method for the neutron capture reac- 
tion 155 Gd+n -> 156 Gd*. Thus, one has to study the 
surrogate reaction for excitation energies E ex of 156 Gd 
that correspond to few MeV neutron absorption in 155 Gd 
(0 < E n < 2); by energy balance one has: E ex = 
Mi55 Gd + M n + E n - Mi56 Gd S n + E n where M is 
the corresponding rest mass and S n = 8.536 MeV is the 
neutron separation energy for 156 Gd pQ. Therefore, we 
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FIG. 1: Reaction cross section for the desired reaction a+A — > 
B* — > C+c might involve an unstable target A. Alternatively, 
the cross section may be deduced using theoretical modeling 
and experimental data from a surrogate reaction d+D — > b + 
B* — > C+c on a stable target D. Both reactions proceed via 
the formation of the same intermediate compound system B*. 

shall focus on excited states in 156 Gd such that: 

E™ Gd > 8.5 MeV. (1) 



II. MATHEMATICAL FRAMEWORK 

Atomic nuclei exhibit a multitude of spectral phenom- 
ena. Among the simplest and relatively well understood 
properties are the rotational and vibrational spectra as 
well as the general single-particle shell structure that ex- 
plains the observed magic numbers [2j [3] . The traditional 
modeling of deformed even-even nuclei considers the ro- 
tational states of nuclei as a collective phenomenon whose 
states are simple rotational functions [3]. An extra par- 
ticle is then assumed to occupy a valence single particle 
state within a deformed mean-field potential. This allows 
one to treat the odd-even deformed nuclei within the ro- 
tor plus single particle model. Within this framework one 
can calculate cross sections for particle transfer as transi- 
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tion amplitudes between collective rotational states and 
a particle plus rotor states. 



A. States of axially deformed nuclei 

1. Bohr-Mottelson rotor model for deformed nuclei 

Within the rotor model an intrinsic state of a deformed 
system, with axial symmetry, gives rise to a rotational 
band with energies and wave functions of the form [3]: 

J = K,K+(l + 6 K ,o),K + 2(l + 6 Ki0 ),... 
E{J) = E K + ^(J(J + 1) - K(K + 1)), (2) 



1 /J(J+1) 

^KJM ~ A \ n x rX 

4tt y (1 + 6k,o) 

x [D J MK (u)$ K + (-I)^X-kH^] • 

Here X is the moment of inertia of the system, which for 
simplicity is assumed to be independent of K, uo stands 
for the angles 6 and </>, D J MK are the rotational matrices, 
&k is the reflection of the intrinsic state <&k with respect 
to the symmetry axis, and K is the angular momentum 
projection onto the symmetry axis. 



B. Single particle states in axially deformed nuclei 

Since the nucleon density in nuclei is constant, except 
near the nuclear surface, one expects that the general 
form of the mean-field potential, which defines the single 
particle states, is of Woods-Saxon type [4]: 



V(r) = 



l + ex p(^) 



(5) 



For spherical nuclei R is constant and represents the po- 
sition of the nuclear surface while a is related to the dif- 
fuseness of the potential near the surface. For deformed 
nuclei, however, R depends on the surface point of inter- 
est and is often parametrized using spherical harmonics: 



R(0, 



Ro 1 



XjJL 



(6) 



For axially deformed nuclei R does not depend on (f) due 
to the axial symmetry. The A = 1 is absent due to cen- 
ter of mass considerations and odd A terms are usually 
absent due to parity considerations; for axially deformed 
nuclei one usually assumes quadrupole and hexapole de- 
formation only: 

R{6, <t>) = Ro(l + f3 2 Y 20 (6, <j>) + f3 A Y 40 (9, <j>)) (7) 



2. Rotor plus particle/hole system 

In the zero order approximation where, we neglect pos- 
sible particle-core couplings such as Coriolis coupling, the 
combined system of a particle/hole plus a core has an in- 
trinsic state that can be viewed as a direct product 
of the intrinsic state of the core and the single par- 
ticle/hole state ?/v 



*n=^**\ E n = E K + e v . 



(3) 



where Q, = K + v due to axial symmetry. Therefore, for 
the particle plus core system we have: 

J = O, n + (l + J n , ), + 2(1 + <y n ,o),» • 
E(J) = Efi + ^(J(J + 1) - + 1)), (4) 



1 / 2J+1 

x [Di fl M§ K i + (-l) n+J D J M _ Q (u;)^ R ^] . 

As long as the final states that we are interested in are 
assumed to be built out of the same initial intrinsic core 
state, we can make perturbation adjustments, presum- 
ably small, to the single particle energies to incorporate 
BCS pairing and Coriolis coupling effects. 



1. Single particle states at small deformation 

For small values of the deformation parameters /3 < 0.3 
in ([7]) one can expand the Woods-Saxon potential ([5| in 
Taylor series. The simplest approach is to consider only 
the first order terms in the expansion [5]. 

V(r, R) ~ V(r, R ) - V'(r, R )^(f3 2 Y 20 (6, <j>) + ...) (8) 

a 

While in many cases this seems to be sufficient, some 
particular single particle states are sensitive to small de- 
formation values of beta fa < 0.1. 



2. Single particle levels at strong deformation 

If (3 is sufficiently small, so that the Taylor expansion 
converges, one can re-express V(r, R) in multiple expan- 
sion in terms of spherical harmonics [6]. An alternative 
way is to solve numerically the Schrodinger equation for 
the deformed Woods-Saxon potential [7]. When compar- 
ing the single particle energies calculated numerically to 
the one calculated using only first order Taylor expansion 
approximation, one finds that for rare-earth nuclei (nu- 
clei near Gd) the "Nilsson diagrams" agree for fa ^ 0.1 
but start to deviate at larger deformations; in particular, 
there is a substantial deviation for rrij = j states. 

To understand the impact of the deformation and what 
should be considered as a small deformation, calculations 
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FIG. 2: 157 Gd deformed single particle energies as function of 
the deformation f3 for F =-44.6883 MeV, r =1.27 fm, a =0.67 
fm, Via = 15.9945 MeV. The circle indicates the region of agree- 
ment between the linear treatment and the exact numerical 
treatment. The arrows at (3 = 0.3 indicate the most sub- 
stantial deviation, between the linear treatment and the ex- 
act numerical treatment, observed for the j = m — 9/2 and 
j = m = 7/2 states. 



were performed using several codes [5j [7J [8] [9j [TQl [TTJ [12] . 
First we made sure that in the case of zero deformation 
the low energy spectrum agrees well within few keV. In 
most cases an initial disagreement was fixed by improv- 
ing the treatment of the reduced mass of the particle-core 
system. We also found that the choice of the spin-orbit 
interaction could result in a factor of about 2 or 4 depend- 
ing on the convention used. Usually obtaining the pi/2 
and p 3 /2 splitting is the main indication that the factor 
affecting the strength of the spin-orbit term is properly 
accounted for. The energy of the Si/ 2 state was used 
for understanding the convention for the strength of the 
central potential. When the energies of the Si/2,Pi/2 5 
and p 3 /2 states are reproduced, by each of the codes and 
the reduced mass is properly accounted for, the rest of 
the bound spectrum is usually well reproduced within 
few keV. For small non-zero deformation we used fewer 



codes [5j [7J [10] , while for large deformation we relied 
only on [7]. 

Fig [2] shows the dependence of the single particle en- 
ergies of the nlj = lhg/2 and 2/7/2 orbitals. Calcula- 
tions using a code that treats the interaction as linear in 
the deformation ([8| are compared to an exact numerical 
treatment of the interaction [5j [7] . Fig [2] illustrates that 
one has to treat deformation with caution. 

In Tabled are shown the neutron bound states in 157 Gd 
calculated using deformed Woods-Saxon potential [15] . 
Since we have used WSBETA code [7 , which uses an ex- 
act numerical method instead of a linear approximation, 
we have set ^4 = which gives the binding energy of the 
47 th neutron single particle state close to the experimen- 
tal neutron separation energy (S n = 6.3598 MeV [14]) in 
157 Gd e 47 = -6.361 MeV. 
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TABLE I: Neutron bound states in 157 Gd calculated with the 
WSBETA code [7 using Woods-Saxon parameters from Ref. 
[T3] but f3 4 = so'that the 47 th neutron state (e 47 = -6.361) 
is near the experimental neutron separation energy S n = 
6.3598 MeV [14], V =-45.1776 MeV, r =1.25 fm, a =0.65 fm, 
Vi s = 19.2015 MeV, p 2 = 0.29 and f3 4 = 0. The Fermi level and 
the level that results in excitation near the neutron separa- 
tion energy in 156 Gd (S n = 8.536 MeV) are in boldface and 
marked by two horizontal lines. 
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C. Adjustments to the single particle energies 



In the Distorted Wave Born Approximation (DWBA) 
the single particle wave function of the transferred par- 
ticle is of essence along with the probability amplitude 
of the particle to be in that state. Ideally one would 
calculate the many body eigenstates of the initial and fi- 
nal system within a chosen single particle basis and then 
would deduce the probability amplitude - the spectro- 
scopic factor for the single particle to be in any of the 
single particle states. Such approach has the advantage 
that a many body Hamiltonian that describes the initial 
and the final systems is used to deduce the single par- 
ticle wave function of the transferred particle. However, 
finding the many body Hamiltonian that works equally 
well across the nuclear chart has proven elusive. 

The absence of clearly successful and unique many 
body Hamiltonian and method for solving the many body 
problem justifies the use of simpler models for deducing 
the single particle wave function of the transferred par- 
ticle such as the particle plus core model described in 
section III Al There are some issues with such models. 
For example, the total anti-symmetrization, pairing in- 
teraction, and particle-core interaction effects are points 
to be addressed. The effect of the anti-symmetrization 
is usually taken into account within the DWBA calcula- 
tion by using the appropriate combinatorial factors [T5] . 
The particle-core interaction is usually dominated by a 
Coriolis coupling and the pairing effects within the core 
are accounted for through the BCS pairing model. Other 
interaction such as coupling to the quadrupole vibrations 
of the core will not be considered here [16]. 



1. BCS pairing 

The pairing interaction is a many body phenomenon 
which is effectively an interaction between pairs of time 
conjugated particles. Exact solutions of the pairing and 
its variations are well known [TTJ [18] [19]. However, an 
approximation know as BCS paring is very convenient 
[20] . Since the pairing interaction is effectively pair cre- 
ation and annihilation process its main effect is on the 
occupation of the single particle states near the Fermi 
energy that represents the location of the last occupied 
state and next unoccupied state. This means that the 
occupation number n v of states that are sufficiently far 
from the Fermi energy \i are practically not affected by 
the paring. 



/i 



V(e,-/i) 2 + A2 



(9) 



Fig. [3] shows the occupation number n v by applying 
the BCS equations to the set of single particle energies 
shown in Table [I] From the graph one can see that the 
occupation of the particle levels with energies e v < —9 
MeV are practicly unaffected by the pairing interaction. 




e v [MeV] 



FIG. 3: Occupation number n v for single particle bound 
states e v in the 157 Gd system. Single particle energies 
e v are calculated for deformed Woods-Saxon potential with 
parameters: Vb =-45. 1776 MeV, ro=1.25 fm, ao=0.65 fm, 
V^ = 19.2015 MeV, deformation /3 2 = 0.29 and the following 
BCS parameters: Fermi energy /x = —6.169 MeV, pairing gap 
A = 1.307 MeV, and pairing strength g = 0.15 MeV. 



Notice that according to ([!]) one has to pull out neutrons 
that are at least 8 MeV below the Fermi level: 



e v < -14.9 MeV. 



(10) 



For this BCS calculation, we used all neutron bound 
states in 157 Gd except the 47 th state which is assumed to 
be Pauli bocked to the pairing interaction due to the odd 
nucleon. The Fermi energy \i was chosen to be between 
the last occupied state and the first unoccupied state in 
the system /i = (e^j + e/^)/2. The values of the pairing 
gap A and the pairing strength g were determined from 
the BCS equations. 



2. Particle core coupling 

It is well know that the Coriolis coupling is an impor- 
tant interaction I2T1 l22l l23l l24l . 



H c = -^{I + 3-+I-U) 



(11) 



The proper consideration results in Coriolis band mixing 
and technical complications similar to the shell model 
wave function type approach. It is known that the first 
order perturbation to the energy is non-zero for Q = \ 
bands only and higher order terms are needed to obtain 
the energy shifts in ft > \ bands [2TJ [22] . 

For this study, however, we consider only one intrinsic 
state of the core (the ground state of the target) coupled 
to various neutron holes in the core. This, therefore, 
eliminates any Coriolis band mixing from our model. 
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There is, however, an important particle/hole core cou- 
pling that we estimate by using first order perturbation 
approximation. This gives us the energy splitting be- 
tween the states Q = \K — v\ and Q = K + v. It is 
clear that according to the simple non-interaction parti- 
cle plus rotor model (|3|, we have two degenerate states 
Vt = \K — v\ and Q = K + v. To estimate the splitting of 
these two states, due to rotation, we look at the expec- 
tation value of the rotational energy of the particle plus 
core system in the intrinsic frame |3]): 

(SI = K±U\ (J C or + Jp) 2 \SI = K±U) = 

= (K\J* or \K) + (v\%\v) + 

+2 (tt = K ± v\ J cor -j p \n = K±v) = 
= {K\J^\K) + {u\T p \u)±2Ku. 

Since the first two terms are the same and represent the 
rotational energy of the core and the particle, one expects 
the energy splitting to be: 

AE(Sl ± = \K±v\)=±c^Kv. (12) 

Here c is a coefficient that relates the moment of inertia 
around the symmetry axis z to the moment of inertia 
around x and y axes. One could estimate its value from 
the mathematical expressions for the moment of inertia 
of a rigid ellipsoid as given in the Appendix section |V A| 
Here, however, we prefer phenomenological adjustment 
that reproduces the experimentally observed splitting be- 
tween = and ft = 3 bands in 156 Gd by using the 
experimental moment of inertia X needed to reproduce 
the Q = ground state band. 

D. Reaction cross sections for deformed nuclei 

1. DWBA direct reactions within the Optical Model 
Potential (OMP) theory 

The Distorted Wave Born Approximation (DWBA) 
considers the initial (final) state of the system to be a 
product of initial \A,a > (final \B,b >) intrinsic states 
of the target, projectile, and distorted wave function x ± 
that depends on the relative coordinates between the tar- 
get and the projectile [TTJ [I2j [15] . The reaction cross sec- 
tion is then calculated using an Optical Model Potential 
(OMP) for x ± an d the appropriate nuclear interaction V 
relating the initial state \A,a> with the final state final 
\B,b >. The transition amplitude then depends on the 
corresponding interaction matrix < Bb\V\Aa >. For one- 
nucleon transfer reactions this matrix element depends 
on the initial and final state of the transferred nucleon, 
the particular interaction that connects these states, and 
the overlap of the other nucleons. For a neutron pickup 
reaction, it reduces to: the neutron wave function i/j in- 
side the target, the strength Dq of a zero range 5 force 
interaction with the projectile, and a spectroscopic factor 



/S 1 / 2 showing how much the target A looks like a remnant 
B plus a nucleon in the state ip (< A\^B >). 

2. Separation energy matching method for proper tail of the 
single spherical states 

For transfer reactions that result in low energy excited 
states the correct asymptotic tail of the wave function, 
which is related to the neutron separation energy, is very 
important in neutron pickup since one expects that the 
last (outer most) nucleon is being transferred. Using 
wave functions that have the desired binding energy by 
adjusting the depth of the Woods-Saxon binding poten- 
tial is one of the simplest and usually very successful ap- 
proximations for calculating the reaction cross section. 
In this approach one keeps the geometric factors of the 
binding potential, such as ro and ao of the potential, fixed 
from systematics but changes the depth of the potential 
until a bound state ip with the desired binding energy 
is found. The approach has even been extended to the 
point where the variable parameter is the strength of a 
surface derivative potential. 



3. Sturmian method for form-factors 

Problems with the separation energy matching using 
only one spherical wave function has been a motivation 
for more accurate description of the relevant form factor 
by inclusion of more than one basis state. The spherical 
harmonic oscillator (SHO) is a traditional basis to ex- 
pand a wave function. However, the tail in a finite SHO 
basis expansion is allays wrong. An interesting alterna- 
tive basis to expand a deformed bound state is to use the 
set of bound states of the spherical part of the interaction 
[13] . However, this expansion does not guarantee correct 
tail either. An expansion would provide the correct tail if 
one of the basis states is a bound state with the same en- 
ergy as the state that is being expanded and all the other 
non zero components correspond to basis states with a 
deeper binding energy. 

In a Sturmian basis all the basis states have the same 
tail as the original state that is being expanded in this 
basis. In order to maintain correct asymptotic tail one 
is now facing a problem where one has to find differ- 
ent wave functions and potential strengths that result in 
the same energy and thus the same wave function tail. 
The Sturmian basis method has been utilized before [25] . 
The method in [25] relies essentially on expressing the 
deformed potential as liner in fa which may not be suf- 
ficient in the case of strong deformation. Unfortunately, 
this approach relies on the linear form of the potential 
with respect to the deformation parameters. 

Another, more successful and physically more appeal- 
ing approach is to consider coupled channel (CC) calcula- 
tions [10 . Unfortunately, this code in its current version 
has only fa deformation. In one of the next section we 
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will discuss briefly our comparison to this code and pos- 
sible further developments. 

To illustrate the role of various parameters involved, we 
now look at the incoherent DWBA reaction cross section 
for a particle transfer from a deformed single particle 
state i\) v that can be expressed in terms of transfer cross 
sections on spherical single particle states <p n ij [15j[25]: 

da(J z K z - J f K f] v)=J2 J2( a » v » c » j ^ da niY ( 13 ) 

lj n 

Here cr^iY are the DWBA cross section from a spherical 
state <j) n ij , c™ l i are the expansion coefficients of the state 
ip u in spherical basis states n /j, v v represents BCS occu- 
pation number (y% = n v /2) of the state and a v is the 
Coriolis band mixing amplitude. The spectroscopic fac- 
tor Sij is often used as shorthand notation for the term 
raised to second power. In our calculations, we actually 
consider the generally more appropriate coherent cross 
section by using super-position of basis states <p n ij with 
amplitudes a v v v c^K 

III. NEUTRON TRANSFER REACTION 
RESULTS 

In this section we present the results of our calculations 
on the neutron transfer reaction 15T Gd( 3 He, 4 He) 156 Gd. 
We first discuss the OMP for the in-going and out-going 
channels, then we show the excited states in 156 Gd ac- 
cording to the particle plus rotor model which in this 
case is actually a hole in the core model. Next we briefly 
discus the model space convergence in the determination 
of the amplitudes. After that the absolute cross sec- 
tions and their convergence are discussed along with a 
smeared cr(E) cross section profile near the neutron sep- 
aration energy. Finally, we turn to a discussion of the 
probability distributions P(J 7r ,£'). 

Our transfer cross sections related to the reaction 
157 Gd( 3 He, 4 He) 156 Gd are for one nucleon removal from a 
deformed single particle state ip ei/ from the 157 Gd system 
viewed as a core in a K = 3/2 ~ state. Presumably, the 
final states of the 156 Gd system are rotational states built 
on the intrinsic state Q just like in the particle plus core 
model ([3]). However, instead particles we work with holes 
in the core \A = 156, tt >= \A = 157, K = 3/2" >. 

A. Elastic cross sections and the optical model 
parameters 

Usually, calculations with distorted waves within the 
DWBA use an Optical Model Potential (OMP): 

U = Vf(x )^iWf(x w )-V so - d ^j^t-s 

r dr 

f(r) = r— 7/ hwT. R = r A 1 / 3 

1 + exp ((r — R)/a) 



Table [II] shows the OMP parameters used in our calcu- 
lation. These parameters are taken from the Reference 
Input Parameter Library (RIPL-2 ) [25] . The Ruther- 
ford normalized elastic cross sections for the ingoing and 
out-going channels are shown in Fig. [4] Notice that the 
absence of interference pattern for th^He cross section 
is due to its relatively low energy. 
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TABLE II: OMP parameters, in fm and MeV units, for the 
in-going 3 He+ 157 Gd [27] and the out-going 4 He+ 156 Gd [28]. 
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FIG. 4: Rutherford normalized elastic cross sections for the 
relevant He+Gd reaction. The smoothness of the 3 He cross 
section is due to the lower value of the associated energy. 



B. Excited states in 156 Gd according to the hole in 
the core model 

In calculating a transfer reaction to an excited state 
in 156 Gd, we first construct an intrinsic state in 156 Gd 
as a hole in the 157 Gd core then we consider rotational 
states built on that intrinsic state. We consider total of 
12 members of the related £1 = \K zb v\ rotational bands 
to include all possible excitations with J < 8. 

\ 156 Gd,n= \K±v\ > = iplj 157 Gd,K = 3/2~ > (14) 
E(J*;Sl=\K±v\) = e -e, + ^(J(J + l)+J ± ) 

Here eo is used to set the ground state energy of 156 Gd 
to zero, and S± is the energy shift of the Q = K + v state 
relative to the Q = \K — v\ state. Thus if we set 5- = 
for Q = \K — v\ then 5+ = 2cKv for Q = K + v states. 
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We always keep z/, and K positive when we use them 
as quantum numbers of the states. 

By fitting the first four excited states of the Q = + 
ground state band, the moment of inertia is fixed to be 
H 2 /21 = 13.59 KeV for the rotational bands in 156 Gd. 



The c = 17.741 coefficient for the energy shift (12) of 
the related Q = 3 + band was chosen to reproduced the 
excitation energy of the 3^ state. The value of c is about 
1.65 times the ratio of the moments of inertia T\\/T± 
for rigid ellipsoid with (3 = 0.29. In Table III is shown 



comparison of the experimentally observed low energy 
states in 156 Gd [14] and excitation energies calculated 
considering a neutron hole coupled to a core. 



TTT 


E exp [MeV] 


E th [MeV] 


J n 


E exp [MeV] 


E th [MeV] 







0.004 


3f 


1.248 


1.255 


2f 


0.089 


0.085 




1.355 


1.363 


4f 


0.288 


0.275 




1.507 


1.499 


6+ 


0.585 


0.574 


6 3 + 


1.644 


1.662 


st 


0.965 


0.982 


n 


1.85 


1.852 



TABLE III: Comparison of the experimentally observed low 
energy states in 156 Gd [14] and excitation energies calculated 
considering a neutron hole v — 3/2 _ at £47 = —6.361 coupled 
to the core system 157 Gd with K = 3/2 ~. The first three 
columns represent the Q = + ground state band; the next 
three columns represent the Q = 3 + band. The theoretical 
ground state is not exactly zero because £47 is not exactly 
equal the neutron separation energy S n in 157 Gd. 

In Fig. [5] are shown the parity (it), angular momentum 
(J), and energy (E) of the states that are considered in 
the calculations of the reaction cross sections <j(J 7r ;£^). 
The absence of positive parity states in the 10 MeV region 
is due to the five positive parity single particle states 
24-28 with energy between -16.81 and -14.66 in Table 
[I] These states will produce a parity asymmetry in the 
P( J 71 ", E) distribution in the energy region of interest for 
the surrogate method. 



C. Model space convergence 

In order to compute the transfer cross section one needs 
the single particle wave function and the correspond- 
ing spectroscopic factor. We already discussed our main 
approximations to the spectroscopic factor: we neglect 
the Coriolis band mixing (a v = 1), the pairing effects 
(v v = 1), and consider single particle wave functions that 
are expressed in a Sturmian basis. The idea is to main- 
tain the same exponential tail for each spherical basis 
state as the tail of the deformed state that we are inter- 
ested in. We have to make expansion in a spherical basis 
in order to be able to use existing reaction codes. 

Here we consider the Sturmian approach by using Stur- 
mian spherical basis states. We calculate the spectra of a 
deformed Woods-Saxon potential using standard bound 
states technique code WSBETA [7]. Then for each state 



s 10— 



1 1 $ 1 1 1 1 * 1 1 r 



: ; : : 





Ji2J 



FIG. 5: Distribution of the calculated 156 Gd excited states. 
The sign of the x-axis encodes the parity of the states and the 
magnitude is equal to the twice of the angular momentum J. 



with energy e we find all the Sturmian spherical basis 
states (zero deformation) <p en ij with nlj labels as for a 
spherical harmonic oscillator up to the N max oscillator 
shell. These basis states are constructed with the reac- 
tion code DWUCK4 [llj. For a fixed e and nlj labels 
the code finds a scaling factor for the original spherical 
potential such that (j) € nlj is a bound state of this new 
potential. This scaling factor is then used to recompute 
the state <p en ij within the WSBETA code in the same 
basis where the deformed state ip u has been computed. 
The expansion amplitudes are then calculated as de- 
scribed in the Appendix Section |VB| and then passed 
back to the reaction code to compute the cross section 
for pickup from the deformed state ^ v using the spherical 
basis states (j) en ij' Since the amplitudes have to be 
added coherently, we have used the coupled channel code 
CHUCK3 [H]. 

In order to test the reliability of our amplitudes and 
the corresponding radial wave functions, we have com- 
pared some of our c l J and (j) e ij against calculations with 
the Rost's code [ID] . The agreement is good but not per- 
fect; discrepancies can be attributed to differences in the 
implementation of the non-zero deformation. 

In order to understand how well a state i(j u is expressed 
in the corresponding Sturmian basis one has to study the 
norm of ty v in the Sturmian basis and how this norm 
converges with the size of the model space. From Fig. [6] 
one can clearly see the model space convergence. 

The state norm converges even faster. Beyond N max = 
10 the smallest non unit norm is only 0.98. Unfortu- 
nately, there are some particularly troubling states. As 
seen from Fig. [7| there are states that have systemat- 
ically smaller norm then the other states and some are 
even getting norm 1.0009 due to numerical issues. 
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FIG. 6: Model space convergence of the Sturmian basis ap- 
proach. The average number of components is the sum of 
the non-zero Sturmian basis components for the 47 occupied 
bound states divided by 47. 
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FIG. 7: Norm convergence as function of the model space. 
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FIG. 8: Maximal deviation between the actually calculated 
eigenenergy of a basis state and the desired energy. Top 
graph: E-deviation for the model spaces 4 < Nm ax < 10. 
Botom graph: E-deviation for the spaces 9 < N max < 16. 



D. Neutron pickup cross sections 

The individual cross sections, for neutron transfer from 
a deformed state ^ v that results in a final state with 
J 77 and energy E as shown in Fig. [5j are calculated as 
coherent cross sections with the code CHUCK3 [12 using 
the amplitudes & nl - times a Clebsch-Gordan coefficient 
and other appropriate factors [25] : 



It is possible that these troublesome features are actu- 
ally result of a deeper numerical problem inherited in any 
standard large basis bound state technique. The prob- 
lems are likely to be related to the fact that only the low- 
est states are well converged and thus possess the correct 
energy; because of that Sturmian states that are higher 
in n are less accurate, which is seen from the deviation 
of their actual energy from the desired Sturmian energy 
shown in Fig. [8] As seen from the top graph in Fig. [8] all 
basis stats computed with WSBETA N max < 10 have de- 
viation within few keV. However, for larger model spaces 
Nmax > 9 there are basis stats whose deviation is getting 
much too big. This could become one of the main issues 
for the presented computational method. 




■f^D x c v nlj x (J f K f \jm, Ji, Ki). (15) 



Here Do is the strength of the zero range transfer poten- 
tial (Dq5(x), Dq=1S [11). The resulting reaction cross 
sections denned for each single particle state a\(e u ^J 7r ) 
are shown in Fig. [9] In what follows the index A will de- 
note the equivalent pair of indexes (e u , J n ). The size of 
the circles reflects the size of the corresponding cross sec- 
tions. As seen from Fig. [9] there is often a good overlap 
of the circles for the model spaces shown. This indicates 
the convergence of the calculated cross sections. 

The cross sections that one can calculate within the 
presented framework correspond to sharp final states. In 
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Excited States in 156 Gd 

a(f,E) for 157 Gd( 3 He, 4 He) 156 Gd 



~5l 



T" 



o 

; °o 
■ : §81 



o 


N m 


ax <=io 




N m 


ax <=ll 




N m 


ax <=12 





N m 


ax <=13 


o 


N m 


ax <=14 




N m 


ax <=15 



• o 



6 



° o 

T a V 



-4 4 

7C2J 




12 16 20 




FIG. 10: Total neutron cross section cr(E) in mb units 
smeared with Lorentzian using energy dependent width F = 
a-\-bE. Dashed line represents F = 0.1 + 0. IE 1 and continuous 
line represents F = 0.01 + 0.0LE. 



FIG. 9: Individual cross sections for neutron transfer from a 
deformed state that results in a final state with J 71 and energy 
E as shown in Fig. [5] The size of the circles reflects the size 
of the corresponding cross sections. 



reality there are widths associated with the single particle 
states as well as with the final states. In order to produce 
a smooth total cross section as a function of the excitation 
energy of the 156 Gd system we consider a Cauchy type 
smearing distribution function (Lorentzian) and define a 
smooth cr(E) cross section [25] : 



1 



4r 



PX ^ 2kA(E-E v ) 2 + F 2 
a(E) = ^ Px (E)cj x , F = a + bE. 



(16) 



In FigjlO] are shown the effects on the total cross sec- 
tion <j(E) of different smearing width choices F. The 
choice of F should eventually be fixed by comparison of 
the theoretical and experimental cross sections. 

For each choice of a model space N max one can com- 
pute a set of cross sections (J\(e Vl J n ). If the model space 
is sufficiently big so that convergence has been archived 
for the c™ l i then one expects that the corresponding cross 
sections would be the same for the converged spaces. It 
was show in Fig. [6] that convergence requires N max > 12. 
Therefore, as expected, the smeared total cross sections 
for N max > 12 are the same as seen in Fig. 11 




FIG. 11: Smeared total cross section for F = 0.1 + 0.1E as 
a function of Nm ax . The identical a(E) for Nm ax > 12 (red 
curve) illustrate the convergence of the results. 



It is clear that if one sums over all possible J n quantum 
numbers then one gets one, that is, J7T P( J 71 "; E) — 1. 



In Fig. 



12 



are shown the formal probabilities P(J n ; E) 
77 state of 156 Gd at energy E through direct 



for exciting , 

neutron pickup via 42 MeV 3 He on 157 Gd target. 

In Fig. 13 are shown P(J n ;E) distributions that are 



E. The P(J*,E) distributions 



of interest to the surrogate method. These graphs show 
that at least in principle there are no dead channels due 
to zero P(J 7r ) and thus any channel that goes through 
J < 8 should be able to provide information on the decay 
probability of the compound system. The apparently 
zero values of the distribution for J > 6 are due to the 
model space used that focuses mostly on J < 8 excited 
states in 156 Gd. 



The smeared cross sections, introduced in the previous 
section, can be used to define the probability to excite a 
state with quantum numbers J 71 ": 

P{J*\E) = ^E^Wa(£>a (17) 



IV. CONCLUDING REMARKS AND 
DISCUSSIONS 

The present computational method is clearly a viable 
approach to particle pick-up reactions in deformed and 
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Probability distrinution at E=8.5 MleT 



FIG. 12: All P(J™;E) distributions for J < 8. 



strongly deformed systems. However, for stripping reac- 
tions the method would need to incorporate the effects 
of the pairing interaction since adding a nucleon would 
involve mostly levels near and above the Fermi level and 
even unbound states and resonances. In heavier nuclei 
the treatment of the valence single particles states above 
the Fermi level would probably face a problem due to the 
treatment of weakly bound states that may render the 
current computational technique useless; this could eas- 
ily be monitored through evaluation of the norm of the 
bound states ip £m and the quality of the Sturmian ba- 
sis (penijm- An alternative to the current approach that 
uses a standard bound state computational technique to 
compute the overlap of the state ip £m and <j) £n ijm would 
be to extend Rost's code to include higher multipoles 
beyond the fa or to solve numerically for i/j ern as gener- 
alized eigenvalue problem for non zero deformation using 
the Sturmian basis states <p £ nijm- 

The surrogate method (Fig. [I]) assumes that the surro- 
gate reaction populates J 77 states of the intermediate nu- 
clear system within the same energy range as the desired 
reaction. The method would certainly fail if the distribu- 
tion P(J n ,E) becomes zero for some relevant values of 
J 71 at the relevant range of energies E since there would 
be no way of deducing the decay probabilities g x (J^ ,E). 
If P(J 7r , J E) = for a given J 71 " at an isolated energy 
E one could perhaps devise an analytic continuation for 
g x (J w ,E). Our calculations, however, show that within 
the assumptions and computational modeling the reac- 
tion 3 He+ 157 Gd 4 He+ 156 Gd* has a well behaved for- 



mation probability P(J*,E) (see Fig. [12] and Fig. [13) 
within a wide energy range relevant to the desired re- 
action 155 Gd+n — > 156 Gd*. Therefore, given an experi- 
mental input on the decay probability P§ x into an exit 
channel \ within the surrogate formation channel 5, one 
should in principle be able to determine g x (J n ,E). 

Experimental input is essential for the fine tuning of 
the model parameters such as single particle Woods- 
Saxon potential, optical model potential for the incoming 
and out going projectiles as well as the choice of T(E) 
used in the smearing function p(E). Above all a compar- 
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FIG. 13: P( J 7 *; E) distributions for energies near the neutron 
separation energy in 156 Gd. From top to bottom the energies 
are 8.5, 9, and 9.5 MeV. The sign of the horizontal coordinate 
corresponds to the parity 7r and its magnitude gives J. 



ison with experiment should be the true measure of the 
applicability the surrogate method. 
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APPENDIX 



can evaluate the overlap matrix of these states: 



A. Q± splitting due to particle-core coupling 

A more rigorous analysis in the intrinsic frame of the 
rotational energy of the hole in the core system gives: 

<n±\j-T- j\n± >= 

= 2#c - ^-(< K\I 2 \K > + < v\j 2 \v >) 

ZX I 



'2J, 



M\ / M\ 



(18) 



Here Q± = |K±z/|, J is the intrinsic angular momentum 
of the combined system (hole and core), / for the core, 
and j for the hole. He is the Coriolis coupling (11). T 



is the moment of inertia of the combined system. T is 
diagonal in the intrinsic frame with T x = T y = 1/(2X±) 
andT, = 1/(2X||). 

For the shape and moment of inertia of rigid ellipsoid 
of mass M one has: 



ii = T E % - R ( 9 > & = R °( x + p y ™(°> *)) 

(19) 



Which gives: 



In 



32tt + /? (25/? + 8^) 
2 (a/5/?- 40F) 2 



(20) 



B. J amplitudes in a Sturmian basis 

As we already mentioned, for strong deformation one 
needs to consider either the full Taylor expansion [6 or to 
use numerical methods to solve the Schrodinger equation 
0. Given the code WSBETA [7] along with the code 
DWUCK [29] plus a few post-processing codes one can 
solve the Sturmian problem and find expansion of any 
deformed state i/j £m ((3 ^ 0) in terms of non-orthogonal 
Sturmian states <j) £n ijm with correct tails given by their 
binding energy e. 



nlj 



(21) 



For this purpose, one solves the following set of eigenvalue 
equations: 



Hf3=o4>£nljm 



(22) 



Since WSBETA gives solutions to these problems in the 
same cylindrical harmonic oscillator basis \NN Z AQ) , one 



NN Z A 

^ = b£ NzA \nn z aq) 



NN Z A 

Here the quantum numbers iV, N Z ,A, and m represent 
the total number of oscillator quanta, the number of os- 
cillator quanta in the z -direction, and the projection of 
the angular momentum and the total angular momentum 
onto the symmetry axis z. Since Sturmian states <fi £n ijm 
with different Ijm quantum numbers are ort honor mal 
due to orthogonality of the spin and angular momen- 
tum components, we have non-orthogonality only within 
a fixed Ijm sub- set: 



Ijm 
^nk 



^snljm 



\0Ekljm) ~ ^2 D nlfm AD t 



NN Z A 
kljm 



NN Z A 



\NN Z A tdNN z A 



l\ z l\ ni' 

/ j nljrn rr 
NN Z A 



(23) 



(24) 



In order to determine the coefficients c nZj in ( |2l| we have 
to use the inverse ^ m of the overlap fi 1 ^™ (J2 k l^nkl^ n ' k = 
5™ ). Within each sub-space of fixed Ij we have: 



E 



C l^nk 



fixed Ij => c n 



/i ^k 

k 



In general one expects to have no problems in inverting 
fi n k since det \i = would imply liner dependence of the 
vectors ^enijm which should not be present since each of 
these states has different number of nodes. How many 
Sturmian states (j) en ij are needed is determined by the 
normalization condition: 



eQ\ 



= ^EE 



= EE c "" s « 



lj nk 



Therefore the number Ns of Sturmian states (j) £ nij is de- 
termined by the normalization requirement: 



J max Ns 

lj n,k 



(25) 



Here an implicit lj index is assumed for E n and fi nk . 
Since parity considerations determine the value of I for 
any given j one can write all lj indexed partitions as 
single j index. 



When the asymptotic tail is proper (22) and Ns and 
jmax are such that the normalization condition (25) is 



satisfied then one can be sure that there are no states 
omitted, continuum or discrete, because we have nor- 
malized bound state with the correct expansion at small 
r as well as at large r. 
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C. Phase conventions for the wave functions 

The overall phase of a wave function, single particle 
wave function in particular, is arbitrary and unobserv- 
able. Thus, it is usually fixed for convenience by a suit- 
able choice of convention. Although there could be as 
many conventions as wave function, there are two main 
choices in the literature. One convention is to consider 
wave functions that are positive as r — > oo [TTJ [12] . An- 
other convention is to consider wave functions that are 
positive at the origin r — » [5] . Sometime the convention 



is not clear at all [7]. 

For our calculations, it is important to adjust the cj^- 
since we use two different codes [TJ [12]. First the c v nl - 
are calculated from the WSBETA 7 wave functions as 



described in Appendix section VB Then we take into ac- 
count the phase difference relative to the wave functions 
generated by DWUCK4 and CHUCK3 [TEJ[12]. For this 
purpose, we evaluate the sign of the WSBETA spheri- 
cal wave functions that are given in cylindrical basis at 
# = 7r/3,(/> = 0, and r — > oo using an asymptotic expres- 
sion of the cylindrical basis wave functions. 
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